Numerical study of the one-dimensional quantum compass model 



o 

(N 



Department of Physics, 



Saeed Mahdavifar 
University of Guilan, 41335-1914, 
(Dated: July 26, 2012) 



Rasht, Iran 



The ground state magnetic phase diagram of the one-dimensional quantum compass model (QCM) 
is studied using the numerical Lanczos method. A detailed numerical analysis of the low energy 
excitation spectrum is presented. The energy gap and the spin-spin correlation functions are cal- 
culated for finite chains. Two kind of the magnetic long-range orders, the Neel and a type of the 
stripe-antiferromagnet, in the ground state phase diagram are identified. Based on the numerical 
analysis, the first and second order quantum phase transitions in the ground state phase diagram 
are identified. 

PACS numbers: 75.10.Jm Quantized spin models;75.10.Pq Spin chain models 



I. INTRODUCTION 
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A theoretical understanding of the magnetic proper- 
ties of low-dimensional quantum magnets has attracted 
a lot of interest in the last decades. Especially, ground 
state properties are interesting since the quantum fluctu- 
ation often plays the dominant role at zero temperature. 
Recently, theoretical works are focused on the important 
role played by the orbital degree of freedom in deter- 
mining the magnetic properties of materials. The com- 
plex interplay among the spin and orbital degrees of free- 
dom in quantum magnets, such as Dzyaloshinskii-Moriya 
(DM) interaction^, makes their phase diagram rich and 
induces various fascinating physical phenomena. 

The experimental observations on Mott insulators are 
a realization of the effect of the orbital degrees of free- 
dom on the low-energy behavior of a system. A good 
candidate for explaining the low-temperature behavior 
of some Mott insulators is the quantum compass model 
(QCM)2. In this model, the orbital degrees of freedom 
are represented by (pseudo) spin- 1/2 operators and cou- 
pled anisotropically in such a way to mimic the com- 
petition between orbital ordering in different directions. 
The two-dimensional QCM is introduced as a realistic 
model to generate protected cubits^ and can play a role 
in the quantum information theory. This model is dual to 
studied model of superconducting arrays^. It was shown 
that the eigenstates are two-fold degenerate and to be 
gapped^. The results from both spin wave study and 
exact diagonalization have suggested a first-order quan- 
tum phase transition in the ground state phase diagram^. 
Recently, the existing of the first-order phase transition 
is confirmed in two-dimensional QCM££. The dilution 
effect is studied by means of the quantum Monte- Carlo 
method^. It is found that due to dilution, the ordering 
temperature decreases much stronger than that in spin 
models. A numerical monte carlo simulation also is done 
on the square lattice and critical temperatures are ob- 
tained for the classical and quantum versions^. 

The one-dimensional (ID) QCM has been studied 
much lessiiriSi. The Hamiltonian of the ID QCM on 
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FIG. 1: (Color online.) The phase diagram of the one- 
dimensional quantum compass model. The first-order and 
the second-order critical lines denote with Ji/Li = and 
J2/L1 = 1, respectively. The short-range correlation func- 



tions in different regions are given by (fji— i ff M )• ^ n ^ ne 
regions III. and IV., there are the stripe AF and Neel long- 
range orders, respectively. 



a periodic chain of N sites is given by-i£ 



N/2 



J 20%, 



2CT 2l _ 1 cr 2i 



2i a 2i+l- 



(1) 



Here &l' 7 ' are the Pauli operators on the zth site and 
Ji,J2,Li are the exchange couplings. By mapping the 
model to a quantum Ising model, an exact solution is 
obtained for the ground state energy and the complete 
excitation spectrum-^. It is shown that the ID QCM ex- 
hibits a first-order phase transition at J\ = between two 
disordered phases with opposite signs of certain local spin 
correlations. The degeneracy of the ground state energy 
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FIG. 2: (Color online.) The Schematic picture of the canted 
spins in the direction of the z axis, which is known as the 
saturated stripe-antiferromagnetic order. 



at this transition point is 2 x 2 N ' 2 in the limit of N 
and the energy gap is closed. The model is also diagonal- 
ized exactly by a direct Jordan- Wigncr transformation-^. 
The obtained results by latter approach, confirm the ex- 
istence of the first-order phase transition in the ground 
state phase diagram. In a very recent paper, Eriksson 
et.al.ji£ have shown that the reported first-order phase 
transition, in fact occurs at a multicritical point where 
a line of the first-order transition {J\/L\ = 0) meets 
with a line of the second-ordcr( J2/L1 = I) transition. 
Also, there are four gapped phases in the regions: (I.) 
J1/L1 < 0, J 2 /Lj > I, (II.) A/Lx > 0, J2/L1 > f, (III.) 
J1/L1 < 0, Ja/Li < I, (IV.) h/Li > 0, J 2 /Li < 1 
(see FigJTJ. Since the spin ordering in different sectors 
is obtained by study of short-range correlation functions 
(( a 2i-i a 2i Z ))i th° complete picture of different long-range 
orders of the ID QCM is unclear up to now. 

For the first time, in this paper we perform an accurate 
numerical experiment on the ground state phase diagram 
of the model. By analyzing the numerical results on finite 
chains, we draw a clear picture of the different long-range 
orders in the ground state phase diagram. In particular 
we apply the Lanczos method to diagonalize numerically 
finite (N = 12, 1 6, 20, 24) chain systems. Using the ex- 
act diagonalization results, we calculate the various spin 
structure factors. Based on the exact diagonalization re- 
sults we argue that instead of the analytical suggested 
Neel order, the stripe-antiferromagnetic long-range or- 
der (see Figj5J) exist in the region (III.) J\jL\ < 0, 
J2/L1 < f of the ground state phase diagram. We de- 
note by "stripe-antiferromagnetic phase" the phase with 
the opposite magnetization along the z axis on the odd 
bonds. 

The paper is organized as follows. In the forthcom- 
ing section we present numerical results of our exact di- 
agonalization studies of the system. In Section III, we 
conclude and summarize our results. 



II. NUMERICAL RESULTS 

In this section, to explore the nature of the spec- 
trum and the quantum phase transition, we used Lanc- 
zos method to diagonalize numerically chains with length 
up to N = 24 and different values of the exchanges (a) 
Ji/Li = 2(1 - Jj/Li), (b) J 2 /Li = f - J1/L1 and (c) 
J2/L1 = (1/2)(1 — Ji/Zi). The energies of the few lowest 




FIG. 3: (Color online.) Difference between the energy of the 
first excited state and the ground state energy as a function 
of the exchange J\, for chain exchanges Li = 1, J2/L1 = 
f — J1/L1, including different lengths N = 12, 16. 



eigenstates were obtained for chains with periodic bound- 
ary conditions. The Lanczos method and the related re- 
cursion mcthod a 17 ' 18 , possibly with appropriate imple- 
mentations, have emerged as one of the most important 
computational procedures, mainly when a few extreme 
eigenvalues are desired. 



A. The Energy Gap 

To show that the transition lines could be observed 
from the numerical calculations of small systems we 
start our consideration with the energy gap function. 
First, we have computed the three lowest energy eigen- 
values of chains with different values of the exchanges 
J2/L1 = 1 — J\jL\. To get the energies of the few lowest 
eigenstates we consider chains with periodic boundary 
conditions. In Fig 131 we present results of these calcu- 
lations for the exchanges L\ = 1, — 2 < J\jL\ < 1 and 
chain sizes N — 12, 16. We define the excitation gap 
as a gap to the first excited state. As it can be seen in 
Fig 131 in the considered limit of exchanges, this difference 
is characterized by the indistinguishable (within the used 
numerical accuracy) dependence on the chain length and 
shows an universal linear decrease with increasing ex- 
change J\jL\. At J\jL\ = —2 the spectrum model is 
gapped. By increasing the exchange J\/ L\ from —2, the 
energy gap decreases linearly with J\jL\. The energy 
gap vanishes at 3\jL\ = 0, which is the only level cross- 
ing between the ground-state energy and the first excited 
state energy. With more increasing J\/L\, the spin gap 
opens again and for a sufficiently large exchange \J\\/L\, 
it increases linearly. Therefore, in complete agreement 
with the analytical results^, we conclude that in this 
case, there are two different gapped phases in the re- 
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FIG. 4: (Color online.) The spin correlation Wf z (n) as a function of the n, for chain length N = 24 and exchanges L\ — 1.0, 
(a)J 2 /ii = 2(1 - Ji/Li), (b)J 2 /Ia = 1 - Ji/£i and (c) J 2 /Li = (1/2)(1 - Ji/Li). 



gions (I.) Ji/Lj < 0, J 2 /Lt > 1, and (IV.) Ji/Li > 0, 
J 2 /£i < 1. 



and the spin structure factor at momentum q defined by 

JV-l 



S aa {q) = J2 Wf a {n)exp{iqn). 



(3) 



B. The ground state phase diagram 

To recognize the different phases induced by the ex- 
changes Ji and J2 in the ground state phase diagram of 
the ID QCM, we implemented the algorithm for finite- 
size chains (N = 12, 16, 20, 24) to calculate the order pa- 
rameters and the various spin correlation functions. 

Because of two type of links, we have computed two 
kind of two-point (short-range) correlation functions de- 
fined as (ct^-i^ 20 ); (°2i*' )(T 2i+i) ■ Where (...) denotes 
the ground state average. In complete agreement with 
the analytical results^, found four regimes with differ- 
ent short-range correlations. In the region J%/L\ < 1, 
antiparallel ordering of spin z component on even bonds 
(( a 2i a 2i+i)) i s dominated with respect to the antiparallel 
ordering of spin x component on odd bonds ((cfi-i '!*)) 
in the region J2/L1 > 1- On the other hand, in the region 
J1/L1 < the ground state is in the s — N/2 subspacc 
with (cfi-icfi) = 1 and in the region Jx/Lx > 0, is in 
the s = subspace with {<J2%-x a 2i) = — 1- Therefore, 
the first-order quantum phase transition at Jx/Lx = 
corresponds to the flip sign of the z component of the 
pair correlation functions of spins on every first bonds 
(( a 2i-i a 2i)) ■ I 11 following we try to draw a picture of the 
magnetic long-range ordering of the system in different 
ground state phases. 

To find the long-range magnetic order of the ground 
state of the system, we start our consideration with the 
spin-spin correlation function defined by 
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It is known that the spin structure factor give us a deep 
insight into the characteristics of the ground stated. In 
Fig|4]we have plotted W* z as a function of n along the 
three paths. In the case of J2/L1 = 2(1 — J\jL\) we ex- 
pect three different magnetic phases in the ground state 
phase diagram, indicated at the regions (I.): J\ jL\ < , 
(II.): < Jx/Lx < 0.5 and (IV.): 0.5 < Jx/Lx < 1.0. In 
the case of J 2 /Lx = |(1 — Jx/Lx), one additional phase 
is expected in the region (III.): —1.0 < Jx/Lx < 0. The 
selected values of the exchange Jx/Lx in the figure cover 
all regions of the ground state phase diagram. It can be 
seen from this figure (FigUJ) that the z component of the 
spins on odd sites are pointed in the same direction with 
the <j\ and others (on even sites) are pointed in opposite 
direction at Jx/L\ = 0.8. This is an indication for the 
Neel ordering in the region (IV.) Jx/Lx > 0, J2/L1 < 1. 
It is surprising that the numerical results in FigH^c) at 
Jx/Lx = —0.4 suggest that a different novel phase can be 
exist in the region (III.) Jx/Lx < 0, J2/L1 < 1. The val- 
ues of the spin-spin correlation function Wf z in FigQJc) 
at Jx/Lx = —0.4, show that the two spins on the sec- 
ond odd bond are pointed in the opposite direction with 
the a\ (or similarly with the two spins on the first odd 
bond) but two spins on the third odd bond are pointed 
in the same direction with the erf. Clearly seen that the 
mentioned behavior continue in whole of the chain sys- 
tem. Wc denote this novel kind of the ordering in the 
ground state magnetic phase diagram of the ID QCM by 
the " stripe-antiferromagnetic phase" . 

Classically, the effect of the negative exchange Jx/L\ < 
is interesting. In the special case of J2/L1 = 0, the 
Hamiltonian reduces to the alternating ZZ Ising model. 
The ground state of the alternating F-AF ZZ Ising model 
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FIG. 5: (Color online.) The spin structure factor S zz (n) as a function of Ji, for chain lengths TV = 12, 16, 20, 24 and exchanges 
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FIG. 6: (Color online.) The correlation function of stripe-antiferromagnetic order \ zz as a function of Ji, for chain lengths 
N = 12, 16,20,24 and exchanges Li = 1.0, (a)J 2 /Li = 2(1 - Ji/Lx), (b)J 2 /Li = 1 - Ji/U and (c) J2/L1 = (1/2)(1 - Ji/Li) 



has a form |*Gs) = I tt-U-tt-U- •••) wltn the long-range 
order canted spins in the direction of the z axis, as illus- 
trated in Figl2] The ordering of this phase is a type of the 
stripc-antifcrromagnctic phase. Therefore, the order pa- 
rameter of the stripe-antiferromagnetic phase is defined 
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For any value of the J\JL\, the Lanczos results lead to the 
staggered magnetization, M z t = J2f=i (~l) Jcr | = ®- 
and M' p = 0, since the ground state is degenerate and in 
a finite system no symmetry breaking happens. However 
the spin-spin correlation function diverge in the ordered 
phase as N — > oo. We computed the spin structure 
factor S zz (ir) and the correlation function of the stripe- 



antiferromagnetic order parameter given by 



N/2-1 

X ZZ = ( E (- i rK--l+^)K--l+2n+^+2n))-(5) 
n=l 



In FigJSJ we have plotted S zz (q = tt)/(N — 1) as a func- 
tion of Ji/Lx for the chain lengths N = 12,16,20,24. 
The Ji/Li-dependency of the spin structure factor, 
S zz (q = 7r), is qualitatively the same as the staggered 
magnetization, M z t . It can be seen that in the region 
(IV.) J\/Li > 0, J2/L1 < 1, the ground state of the sys- 
tem is in the Neel phase. At J\jL\ = 1.0 (or J2/L1 = 0), 
the saturated Neel phase is dominated in well agree- 
ment with the zz Ising model. Induced quantum fluc- 
tuations by increasing J2/L1 from zero (or decreasing 
J\/L\ from the value one), decreases the staggered mag- 
netization from saturation value. In all cases it is clearly 
seen that the staggered magnetization as a function of 
the Ji/ii, displays a jump to zero at the critical value 
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FIG. 7: (Color online.) The spin structure factor S zz (tt) for 
L\ — 1.0 and J2/L1 = 2(1 — J\/L\) plotted as a function of the 
inverse chain length 1/N for different values of the exchange 
Ji. The dashed line which correspond to Ji = 0.5 and J2 = 
1.0 marks the second order quantum phase transition. 




FIG. 8: (Color online.) The correlation function of the stripe- 
antiferromagnetic order parameter x zz f° r Li = 1.0 and 
J2/L1 — 1/2(1 — J1/L1) plotted as a function of the inverse 
chain length 1/N for different values of the exchange Ji < 0. 
The diminishing behavior of the numerical results (dashed 
lines) in the region (I.): J1/L1 < — 1, J2/L1 > 1, shows that 
the stripe-antiferromagnetic order in the region (I.) is not true 
long-range (local ordering). 



J1/L1 = which is known as the first-order phase tran- 
sition. There is not the Neel ordering along z axis in the 
region J\/L\ < 0. Overlapping of the numerical results 
in all regions, except the region < J\/L\ < 0.5 and 
J2/L1 > 1 (FigJSJa)), shows a divergent behavior of the 
function S zz (ir) by increasing the size of chain N. This 
justifies that the Neel ordering along the z axis is true 



FIG. 9: (Color online.) The spin-spin correlation Wf x (n) as 
a function of the n, for chain length N = 24 and exchanges 
Li = 1.0 and path J2/L1 = 2(1 - J1/L1). 



long-range order in the region (IV.) of the ground state 
phase diagram. To check the existence of the Neel order 
in the thermodynamic limit N — > 00 of the system in 
the region (II.): < J\jL\ < 0.5, J2/L1 > 1 , we have 
plotted in FigG]thc N dependence of S zz (n)/(N - 1) for 
different values of J\/L\ > 0. As is seen from this fig- 
ure in the region < J\jL\ < 0.5, J2/L1 > 1 there is 
a diminishing behavior which shows that the Neel order 
in the region (II.) is local (see also FigHJa)). Therefore, 
we conclude that the Neel order, only exist in the region 
(IV.) of the ground state phase diagram and a second- 
order quantum phase transition happens at critical value 
J2/L1 = 1.0. 

An additional insight into the nature of different phases 
can be obtained by studying the correlation function \ zz - 
In FigEl we have plotted \ zz as a function of the ex- 
change J\/L\ for different values of the chain length 
N = 12, 16, 20, 24. In all cases it is clearly seen that 
there is no stripe-antiferromagnetic ordering along the z 
axis in the region J\/L\ > (regions II. and IV.). In 
contrast, in the region (III.) —1 < Ji/L\ < 0, J2/L1 < 1 
(see Figj6^c)), a profound stripe-antiferromagnetic order 
exist in the z direction. In the limit J\jL\ — > 0~ (or 
J2/L1 — > 0.5 + ), by checking the numerical value of the 
X zz we found that the saturated stripe-antiferromagnetic 
phase dominate in the ground state phase diagram. The 
induced quantum fluctuations by increasing J2/L1 from 
0.5 (or decreasing J\/L\ from zero), decreases the stripe- 
antiferromagnetic order from saturation value. Also, it is 
clearly seen that the stripe-antiferromagnetic order dis- 
plays a jump to zero at the critical value J\ /L± = which 
is also an indication of the first-order phase transition. 
To check the existence of the stripe-antiferromagnetic 
order in the thermodynamic limit N — > 00 of the sys- 
tem in the region (III.): J\jL\ < 0, J2/L1 < 1, we 
have plotted in FigJSJthe N dependence of \ zz I {N/2 — 1) 
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for different values of J\jL\ < 0. As is seen from this 
figure in the region — f.O < J x /L x < 0, J2/L1 < 1 
there is a divergent behavior which shows that the stripe- 
antiferromagnetic order in the region (III.) is true long- 
range order. On the other hand the diminishing behavior 
of the numerical results in the region (I.): J\jL\ < — I, 
J2/L1 > 1. shows that the stripc-antiferromagnetic or- 
der in the region (I.) is local. Therefor, by investigating 
the N dependence of \ zz /(N/2 — 1) for different values 
of J x /L x , we found that there is no long-range stripe- 
antiferromagnetic order in regions (I. and II. and IV.) 
and only long-range stripe-antiferromagnetic order ex- 
ist in the region (III.)- However, the numerical results 
on finite chains plotted in FigEJb) , may suggest that 
there is the stripe-antiferromagnetic long-range order in 
a very narrow region very close to the multicritical point 
J\/ L\ = 0, J2/L1 = 1.0, and at this point jumps from a 
non-saturate value to zero. 

To complete our numerical study of the ground state 
magnetic phase diagram of the model, we have com- 
puted the x component of the spin-spin correlation func- 
tion Wf x in all cases. As an example, in FiglH] we 
have plotted W x x as a function of n along the path 
J2/L1 = 2.0(1 — J1/L1). Clearly seen that along the x 
axis, there are not any conventional magnetic long-range 
order. We also found the same results for other paths. 



ID quantum compass model. We have implemented the 
Lanczos method to numerically diagonalize finite chains. 
Using the exact diagonalization results, first we have cal- 
culated the energy gap and various short-rang correlation 

functions ((&2i- X a 2i Z ' ') an d ( a 2i^ a M+i)) ■ I n com pl e te 
agreement with analytical results^, we have showed that 
there are four regimes with different short-range correla- 
tions, which are separated by lines of the first and second- 
order phase transitions (see FigJT]) . Then to answers a 
very important question: "What are the long-range or- 
dered phases in the ground state phase diagram of the 
system?" , by plotting the spin-spin correlation functions, 
we have found that two kind of magnetic long-range or- 
ders exist in the ground state phase diagram, a type of 
the stripe-antiferromagnetic and the Neel orders in the 
regions (III.) J x /L x < 0, J 2 /L x < 1 and (IV.) J x /L x > 0, 
Ji/L x < 1, respectively. By a detailed analysis of the 
numerical results on the spin structure factors and the 
correlation function of the stripc-antiferromagnetic order 
parameter, we have shown that the Neel and the stripe- 
antifcrromagnet orderings are true lang-range orders. 
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